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ABSTRACT 

We report on the first determination of the distance to the Coma Cluster based 
on surface brightness fluctuation (SBF) measurements obtained from Hubble Space 
Telescope WFPC2 observations of the bright EO galaxy NGC 4881 in the Coma 
Cluster and ground-based observations of the 'standard' El galaxy NGC 3379 in the 
Leo-I group. Relative distances based on the I-band fluctuation magnitude, fnj, are 
strongly dependent on metallicity and age of the stellar population. However, the 
radial changes in the stellar populations of the two giant ellipticals, NGC 3379 and 
NGC 4881, are well described by published Mg2 gradients, and the ground-based 
measurements of mj at several radial points in NGC 3379 are used to calibrate mj in 
terms of the Mg2 index. The distance to NGC 3379, assumed to be identical to the 
average SBF distance of the Leo-I group, is combined with the new SBF measurements 
of NGC 4881 to obtain a Coma Cluster distance of 102 ± 14 Mpc. Combining this 
distance with the cosmic recession velocity of Coma (7186 ± 428 km s" 1 ), we find the 
Hubble constant to be Hq = 71 ± 11 km s _1 Mpc -1 . 

Subject headings: distance scale — galaxies: distances and redshifts — galaxies: 
individual (NGC 4881, NGC 3379) 



1 Based on observations with the NASA/ESA Hubble Space Telescope and partially supported by NASA through 
grant HF-1066.01-94A awarded by the Space Telescope Science Institute, which is operated by AURA, Inc., under 
NASA Contract No. NAS5-26555. 
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1. Introduction 



To determine Hq, one needs to measure precise distances to galaxies sufficiently distant for 
local departures from a smooth Hubble flow to be of negligible importance. Consequently we need 
to extend the distance scale to distances of order 100 Mpc or more. Although the Hubble Space 
Telescope (HST) has greatly extended the range of direct Cepheid-based distances, it cannot 
reach the region where the Hubble flow dominates, and despite the successful determination of 
many new Cepheid distances with HST, the dispute over the global value of the Hubble constant 
is still not settled. Methods to reach the Hubble flow via the Virgo Cluster through the use of 
various secondary distance indicators have tended to obtain values larger than 70 km s _1 Mpc -1 
(Freedman et al. 1994, Mould et al. 1995, Kennicutt et al. 1995| ), whereas methods using a Cepheid 
calibration of the peak brightness of type la supernovae consistently have given values below 60 
km s _1 Mpc -1 ( [5aha et al. 1995 , Iammann & Sandage 1995, Sandage et al. 1996, 3aha et al. 



1996, Schaefer 1996a| , [Schaefer 1996b| , Branch et al. 199~6| ). Recently, methods based on correlations 



between light curve shapes and the luminosity of type la supernovae (|Riess et al. 1995j , [Riess et al 



1996| , [Hamuy et al. 1996|) , as well as the nebular SN la method based on modeling of the forbidden 
line emission at late phases ( Ruiz-Lap uente 1996| ), have resulted in values in the range 60-70 km 
s" 1 Mpc- 1 . 



Part of the problem is evidently the significant depth of the Virgo Cluster ( [Yasuda fc Okamura 
1996 ). This uncertainty is eliminated if secondary indicators can be applied to extend distance 



measurements directly to more remote clusters. Hjorth &: Tanvir (1997) have combined the 
Cepheid distance to the Leo-I group ( [Tanvir et al. 1995| ) with a comparison of the fundamental 
plane (FP) of early-type galaxies in the Leo-I group with the FP of the Coma Cluster (Jorgensen 
et al. 1996| ). Also, Baum et al. (1995, 1997) have directly compared the globular cluster luminosity 
functions for two giant ellipticals in the Coma Cluster with the combined luminosity function for 
the Milky Way and M31. These studies also find values in the range 60-70 km s _1 Mpc" 1 . In the 
same spirit we shall, for the first time, derive the distance to the Coma Cluster based on surface 
brightness fluctuations. 



The SBF Method 



The SBF method was introduced by Tonry & Schneider (1988) as a precise secondary distance 
indicator for early-type galaxies. Until now applications of SBF distances have been limited to 
ground-based /-band ( Tonry et al. 1989 , Tonry et al. 199(i| , Tonry 1991 , Jacoby et al. 1992| , |Tonry 
et al. 1997|) and K-band ( [Luppino fc Tonry 1993| , [Paine fc Mould 1994| , pensen et al. 1996| ) 
observations. The superb resolution of the HST should, however, make the WFPC2 camera an 
ideal instrument for obtaining SBF distances to early-type galaxies in the Coma Cluster, the main 
difficulty being the relatively modest aperture of the telescope. A typical luminosity fluctuation 
must produce at least one detected photo electron during the total exposure time in order not 
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to be overwhelmed by the photon noise. The present work is based on eight 900 s exposures of 
NGC 4881 obtained with the PCI camera in the wide-band F814W (I) filter. The basic data 
reductions are fully described by Baum et al. (1995). Our final SBF flux / is equivalent to ~ 0.7 
detected photons in the i ? 814VF-band in a total exposure time of 7200 s. 

Surface brightness fluctuations are due to the random distribution of discrete point sources 
(stars) across the field of the detector. / is defined as the flux of a hypothetical 'fluctuation star' 
which, when randomly spread across the detector, produces the same flux mean and variance per 
pixel as the observed stellar population. As the statistical properties of the power spectrum are 
completely determined by the mean and variance of the observed flux per pixel, we are allowed 
to simulate luminosity fluctuations by a spatially random distribution of a suitable number (in 
this case ~ 1600 stars/pixel) of 'fluctuation stars' over the field of the detector. We have thus 
obtained an accurate estimate of the random error, as well as the statistical bias, by reducing 
many such simulated fluctuation fields. A high resolution Poissonian stellar fluctuation field is first 



convolved with a point spread function (PSF) created by the Tiny Tim (V4.0) program (Krist 



1993). Subsequently, this array is rebinned to the resolution of the PCI camera, resulting in a 
noiseless stellar fluctuation frame. Finally, Poissonian photon noise and Gaussian read noise are 
added. The whole sequence is repeated 50 times, and all these frames are reduced in the same way 
as the real data. The sample average of the derived SBF magnitudes is 0.05 ± 0.03 mag fainter 
than the true value (i.e., we must subtract 0.05 mag from the derived values), and the standard 
deviation on a measurement of a single simulation (8 x 900 s) is 0.22 mag. A Kolmogorov-Smirnov 



test ( Press et al. 1992 ) did not detect any deviation from normality. Our error estimate relies 
completely on the Monte Carlo simulations, which assume that the PSF is well described by the 
Tiny Tim model. Fortunately, Figure 1 shows that the power spectrum of the Tiny Tim PSF is 
a nearly perfect match to a high signal-to-noise power spectrum of an annular region near the 
center of M32. The power spectrum is obtained from a 400 s F&1AW exposure with the PCI 
camera (HST archive data, ID-5464), and the fluctuation amplitude is derived from the same 
angular area and by the same procedure as we use for the NGC 4881 data. This shows that it is 
fully acceptable to use the Tiny Tim model instead of an empirical PSF. 

We have chosen to derive the fluctuation magnitude mi for an annular region between 2 / .'55 
and 5'.'46 from the nucleus of NGC 4881. In that region, the white noise will be dominated by 
photon statistics, not by the read noise (R ^ 5'.'5), and a good fit to a model galaxy is obtained 
{R ^ 2'.'5). This letter is not the place for a detailed description of the SBF method, especially 
as many of the reduction procedures are quite similar to the ones described by Sodemann & 
Thomsen (1995). Instead, we shall concentrate our description on the steps that are specific for 
this particular application. 

The SBF amplitude, Po = /i is obtained from the radial power spectrum of the normalized 
residual image and the radial power spectrum of a high resolution PSF (a factor of 4 sub-sampling) 
as calculated by the Tiny Tim (V4.0) program ( |Krist 1993; ). The volume of the PSF is normalized 
to one, so the derived SBF amplitude (flux) corresponds to a total magnitude. Specifically, Pq 
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is calculated using equation (15) in Sodemann & Thomsen (1995). We must emphasize that the 
lowest wave numbers are omitted, as they are influenced by flat-fielding and galaxy fitting errors. 

In Figure 2 we plot the scaled power spectrum of the PSF on top of the radial power spectrum 
for the annular region in NGC 4881. The constant power spectrum of the white noise, 0.030 DN, 
has been subtracted from the plotted data points. Note that the high power values with wave 
numbers less than 12 are not used in our derivation of Po- We wish to thank the anonymous 
referee for pointing out that the power spectrum of the NGC 4881 data appeared more noisy 
than the power spectrum of the simulated data, i.e., the power spectrum of the white noise in the 
NGC 4881 data was significantly higher than expected from photon noise and assumed read noise 
(5.2e~). Several identifiable sources evidently contributed to this excess white noise: the actual 
read noise (6.2e~) was higher than had been assumed, there was a Poissonian noise associated 
with the subtracted dark current (5.1e~), and the effective number of frames was slightly less than 
8 due to the rejection of pixels containing cosmic ray events. However, these identifiable sources 
did not completely explain all of the excess white noise. Some of it must have been introduced by 
the pipeline processing. 

In order to evaluate the importance of the excess white noise, we added the missing amount of 
Gaussian noise to 25 of the already existing simulated fluctuation frames. The standard deviation 
on a measurement of a single simulation increased from 0.22 mag to 0.28 mag. Figure 3 shows a 
plot of the scaled power spectrum of the PSF on top of the radial power spectrum for one of the 
simulated fluctuation fields with the correct amount of Gaussian noise added. Although greater 
noise reduces the precision of our SBF distance to the Coma Cluster, it does not alter the SBF 
distance itself. 

Our calculated amplitude is / = 0.0140 DN in 900 s, which according to equation (7) of 
Holtzman et al. (1995) gives a total fluctuation magnitude of m,FsiAw(oo) = 33.61. We must 
add 0.10 mag to this value in order to bring it onto the instrumental aperture (C/5) system 
defined by Holtzman et al. (1995). In addition we must add 0.05 mag because at the time of 
the NGC 4881 observations the CCD was operated at —76 °C instead of —88 °C, which was 
used after 23 April 1994 ( Holtzman et al. 1995| ). We can transform the instrumental magnitude 



^fsuw t° Cousins mi by applying Holtzman et a/.'s (1995) synthetic transformation, which is 
equivalent to (fnj — niF8i4w) = —0.05 ( [Holtzman et al. 1995| , Fig. 9) for a typical red fluctuation 
color of (my — mi) ~ 2.6. Although it is possible that the transformation from the ^8141^ 
band to the Cousins /-band may be different for the fluctuations than it is for stars, we derived 
"Ji?8i4iy(oo) = 22.90 for a similar region in M32, which is consistent with the transformation in 
light of the mi values given for M32 by Tonry et al. (1990) and Sodemann &; Thomsen (1996). 
The fact that the friFsiAW f° r M32 (where (jny — mi) ~ 2.3), transformed according to Holtzman 
et al., agrees with the ground-based mi for this galaxy does not completely ensure that there 
will be similar good agreement for the redder stellar populations of giant ellipticals. In order to 
settle this question we have compared (m; — mp^iAw)-, as derived from stellar population models 



( Worthey 1994 ) having solar metallicity and ages between 2 Gyr and 17 Gyr, with (I — ^81414^) 
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for the giant stars used in the models. The fluctuation magnitude tTifsiaw transforms, within 
a few hundreds of a magnitude, in the same way as the F814W magnitude of giant stars if one 
chooses the stars red enough to match the SBF colors. Actually, a 12 Gyr model, which matches 
the average Mg 2 index of our annulus in NGC 4881, has an offset of (mj — niFmiw) = —0.16 and 
a fluctuation color of (my — mi) = 2.75. A red giant library-star with (V — I) = 2.79 has an offset 
of (J - F8UW) = -0.16. The difference between this value and the offset (7 - F8UW) = -0.05 
given by Holtzman et al. (1995) is almost certainly due to differences in the / passbands used 
to derive the synthetic colors. Holtzman et al. used Landolt UBVRI response curves, whereas 
Worthey (1994) applied the UBVRI curves of Bessell (1990). The difference between Landolt 
& Bessell filter passbands (~ 0.1 mag) is an indication of the remaining uncertainty in the 
transformation of F814W for M stars. In the following we shall use (mi — mpsuw) = —0.10. 

The variance per pixel due to unmasked globular clusters (GCs) can be written as 
°GC = ^GC^Gc/gc ( Sodemann &; Thomsen 1995| ), where ncc is the number of detected GCs per 



pixel, /gc is a fl ux defining the location of the apparent GC luminosity function (GCLF), and Qgc 
is a factor that is calculated from the apparent GCLF and the completeness function (a detection 
probability). For a Gaussian LF, <5gc is given by equation (4) of Sodemann & Thomsen (1995). 
Following the arguments of Baum et al. (1995, 1997) we shall use a model based on transporting 
the MW+M31 GCLF to a Coma distance of (m - M) = 35.2 ( [Baum et al. 1995| ). The combined 



GC data for the Milky Way and M31 are well fitted by a hyperbolic GCLF of the form (Baum 



et al. 1997 ) logn/n = — (a/b)^b 2 + (m — m-Gc) 2 + k(m — mcc), where w-gc is the magnitude of 
the flux, /gc> at the intersection of the asymptotes. As the detection probabilities were derived 
for the V-band the completeness function had to be shifted by the mean color (V — I) = 0.97 



(Baum et aZ.~T995|) of the GC system in NGC 4881. 



In a similar way, the variance per pixel due to unmasked background galaxies can be expressed 
as <7g = Qgn g f GC , where n g is the number of detected background galaxies per pixel and Q g is a 
factor that is calculated from the apparent galaxy LF and the completeness function. The slope 
of the galaxy LF is taken as the average of the values given for the /-band by Tyson (1988) and 
Lilly et al. (1991), respectively. For galaxies it is not so obvious that we can use our shifted 
completeness function for the U-band to calculate Q g for the /-band, because galaxies have 
a much larger color spread than GCs. Fortunately, the relative contribution from galaxies, as 
compared to the one from GCs, is only a g /a GC = n g Q g / (tiqcQqc) 0.05 so we can safely neglect 
the error introduced by assuming that background galaxies have the same colors as the GCs. 

The contribution to the SBF flux from unmasked GCs and background galaxies is calculated 
as A/ = {QGC n GC + Qgn g )fQc/g, where g is the mean galaxy flux per pixel. There are 13 detected 
GCs inside the annulus [3", 5"], and we expect only about 0.5 background galaxies brighter than 
detection threshold within the same annulus ( Baum et al. 199~5| ). Using these values we find that 



we must decrease the measured SBF flux by A/ = (2.26 ± 0.60) 10 -3 DN. We must accordingly 
add 0.19 ± 0.06 mag in order to correct for the undetected point sources. If we instead use the 
Gaussian GCLF (Vgc = 27.8, a = 1.4) suggested by Whitmore (1996) we obtain a correction 
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of 0.17 ± 0.05 mag, which is not significantly different from the value derived for the hyperbolic 
GCLF. 

The K-correction for the integrated spectrum of an elliptical galaxy at the Coma distance 
is quite modest for the /-band, but fluctuation colors are very red, so it is much larger for mj. 
The K-correction for Mfsuw is K{mF%\4w) = (7.1 ± 0.3)z = 0.16 ± 0.01 mag, where the redshift 
(z = 0.0224) of NGC 4881 relative to the Sun ( |de Vaucouleurs et al. 1991] ) has been used. The 



quoted error was estimated from a comparison of the K-corrections derived from stellar population 
models ( Worthey 1994| ) having solar metallicity and ages between 5 Gyr and 12 Gyr. The Galactic 



absorption in the direction of NGC 4881 is Aj = 0.02 (|Burstein fc Heiles 1984| ). The final corrected 
fluctuation magnitude for NGC 4881 is then m/(4881) = 33.62 ± 0.29. 



3. The Distance to NGC 4881 

We shall calibrate the SBF distance to NGC 4881 in terms of the distance ^3379 = 10.7±0.3 Mpc 
(random errors) to NGC 3379 = M 105 in the Leo-I group. We have adopted the average SBF 
distance for the Leo-I group, as recently derived by Tonry et al. (1997) based on the Cepheid 
distance scale. Sakai et al. (1997) obtained a direct determination of the distance (^3379 = 11.5 
Mpc ±0.7 (random errors) ±1.2 (systematic errors)) to NGC 3379 from the tip of the red-giant 
branch (TRGB), but the systematic errors are relatively large. The spiral galaxy M 96 interacts 
with a giant 200 kpc diameter Hi ring (Schneider 1989) orbiting the two central galaxies in the 



Leo-I group, the El galaxy NGC 3379 and the SO galaxy NGC 3384. Thus the distances to 
NGC 3379 and M 96 are tied together by a physical connection. Tanvir et al. (1995) derived a 
Cepheid distance (dM96 = H.6 ± 0.9 Mpc) to M 96 which is very similar to the TRGB distance 
obtained for NGC 3379. A significantly shorter Cepheid distance (d,M95 = 10.0 ± 0.3 Mpc) has 



been derived for M 95 ( Graham et al. 1997 ), but it is not yet clear whether M 95 is at the same 
distance as NGC 3379. 

The prototype of a normal elliptical galaxy, NGC 3379, has both a well determined radial 



Mg 2 gradient Mg 2 = (0.238 ± 0.002) - (0.065 ± 0.003) log(R/R e ) where R e = 56" (pavies et al 



1993), and reliable measurements of the I-band fluctuation magnitude fnj at four radial distances 
between 14'.' 7 and 48'/ 5 flSodemann &: Thomsen 1995| , Table 5). This is important because mi 
is very sensitive to radial differences in metallicity and/or age of the stellar population (Tonry 



1991, Ajhar & Tonry 1994, Sodemann fc Thomsen 1995| , [Tonry et al. 199"7 ) as indicated by 



differences in, e.g., the Mg 2 index. As shown by a plot of mi against (B — I) ( jSodemanri 



fc Thomsen 1995 , Fig. 7) the innermost point, being the most difficult to measure, deviates 
significantly from the expected linear relation. Using a linear regression on mi as a function of 
Mg 2 , excluding the innermost point, we find that m/ = (28.77 ± 0.03) + (7.9 ± 0.7) (Mg 2 - 0.26) 
for the stellar population in NGC 3379. The assigned error on the constant term mainly reflects 
the photometric zero point error. This formula should only be used to interpolate mi in the 
metallicity range 0.24 ^ Mg 2 £ 0.28. 
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Fortunately, Thomsen & Baum (1987) have measured the magnesium gradient for NGC 4881 to 
be AMg 2 /Ai? 1 / 4 = -0.095 ±0.027 mag arcsec ' . Combining this gradient with Dressier 's (1984) 
measurement Mg 2 = 0.292±0.007 of the Mg 2 index within a 4" x 4" region centered on the nucleus of 
NGC 4881 we obtain the following relation: Mg 2 = (0.292±0.007) - (0.095±0.027)(( J R 1 / 4 ) - (i?o /4 )), 
where the mean of i? 1//4 for a disk with the same area as Dressler's slit is (R 1 ^) = 1.089, and 
the mean of i? 1 / 4 for the annulus between 2 / /55 and 5'/46 is (i? 1 / 4 ) = 1.424. Our best estimate 
of the Mg 2 index for the SBF annulus in NGC 4881 is thus Mg 2 = 0.260 ± 0.011. Using our 
relation between rnj and Mg 2 , a K-correction of Kirni) = 0.02, and a Galactic extinction of 
A T = 0.02 ( [Bur stein fe Heiles 1981 ), we find that a stellar population in NGC 3379 with the 



same Mg 2 index as the one in our SBF annulus in NGC 4881 has a fluctuation magnitude of 
to/(3379) = 28.73 ± 0.09. 

Based on our SBF measurements the best estimate of the relative distance modulus between 
NGC 3379 and NGC 4881 is A(m - M) = 4.89 ± 0.30 (random errors only), corresponding 
to a (luminosity) distance ratio of ^4881/^3379 = 9.5 ± 1.3. Using our adopted distance to 
NGC 3379, ^3379 = 10.7 ± 0.3 Mpc, we finally obtain a distance modulus for the Coma Cluster 
of (m — M)o = 35.04 ± 0.31 (random errors only), corresponding to a luminosity distance of 
d 4881 = 102 ± 14 Mpc. 



4. The Hubble Constant 

At small redshifts the distance modulus is given by (to — M)q = 25 + 51og(oz/i?o) + Am(«), 
where Am(z) = 1.086(1 — qo)z QHeckmann 1942] ), z = 0.024 is the redshift of the Coma Cluster 
with respect to the cosmic background radiation ( |Han fc Mould 1992|) , and Hq is the Hubble 
constant in km s _1 Mpc . If we assume that qo = 0.25 then Am(z) = 0.02. We adopt 
cz = 7186 ± 428 km s™ 1 ( [Han fc Mould 1992| ) as likely values for the recession and peculiar 
velocities of the Coma Cluster. Finally, using these values and our distance modulus for the Coma 
Cluster, (to — M)q = 35.04 ± 0.31, we estimate the Hubble constant to be Hq = 71 ± 11 km s~ 4 
Mpc" 1 . 

We are grateful to Jens Hjorth for communication of unpublished results, and for many 
useful comments and discussions. We wish to thank Marianne Sodemann for carefully reading the 
manuscript. 
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Fig. 1. — The scaled power spectrum of the Tiny Tim PSF (line) is compared with the radial 
power spectrum of the surface brightness fluctuations in an annular region near the center of M32 
(squares). The SBF amplitude Pq is derived from wave numbers between 12 and 36 (1.0 — 3.1 
arcsec -1 ), which is the range used for the NGC 4881 data. 
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Fig. 2. — The scaled power spectrum of the PSF (line) is compared with the radial power spectrum 
of the surface brightness fluctuations in the annular region [2755, 5'.'46] centered on the nucleus of 
NGC 4881 (squares). The constant power spectrum of the white noise -Pi (0.030 DN) has been 
subtracted from the plotted data points. The SBF amplitude Po is derived from wave numbers 
between 12 and 36. 



-12- 



<: 
O 



0.06 ~ 
0.04 - 



- □ 

F- □ 



-□ 

nr □ 



0.02 
0.00 
0.02 

0.04 h 

0.06 



"i 1 r 



□ □ 



J I L 







"i 1 r 



□ 



□ 



□ 



J I L 



"i 1 r 



□ 



□ 



J I L 



"i 1 r 



□ 



J I L 



60 



.0 



ave number 



Fig. 3. — The scaled power spectrum of the PSF (line) is compared with the radial power spectrum 
of a simulated fluctuation field covering the same area as used in NGC 4881 (squares). The assumed 
SBF amplitude / has been reduced by the contribution from undetected point sources. The constant 
power spectrum of the white noise has been subtracted from the plotted data points. The SBF 
amplitude is derived from wave numbers between 12 and 36. 



